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Abstract 

Recent investigations suggest that the neutrino-heated hot bubble 
between the nascent neutron star and the overlying stellar mantle of a 
type-II supernova may be the site of the r-process. In the preceding 
a-process building up the elements to A ~ 100, the 4 He(2n,7) 6 He- 
and 6 He(a,n) 9 Be-reactions bridging the instability gap at A = 5 
and ^4 = 8 could be of relevance. We suggest a mechanism for 
4 He(2n,7) 6 He and calculate the reaction rate within the a+n+n ap- 
proach. The value obtained is about a factor 1.6 smaller than the 
one obtained recently in the simpler direct-capture model, but is 
at least three order of magnitude enhanced compared to the pre- 
viously adopted value. Our calculation confirms the result of the 
direct-capture calculation that under representative conditions in the 
a-process the reaction path proceeding through 6 He is negligible com- 
pared to 4 He(an,7) 9 Be. 
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1 Motivation 



Many calculations of the two-step radiative capture reactions of astrophysical 
interest have been done by now in the framework of the single-particle direct- 
capture model. An accuracy of the model was not known to a sufficient degree 
and its assessment by performing a few-body type calculation is of interest. 
This is accomplished in the present work for the 4 He(2n, 7) 6 He reaction. 

In the past years a new neutron-rich astrophysical scenario taking place 
in the neutrino-heated hot bubble between the nascent neutron star and 
the overlying stellar mantle of a type-II supernova has been proposed |Q- 
M. The material, originally in nuclear statistical equilibrium (NSE) at high 
temperature, is expanded and cooled so rapidly that not all the a-particles 
have time to reassemble. In this environment most nucleons are either in the 
form of free neutrons or bound in a-particles. The nucleosynthesis in the 
neutrino bubble takes place in subsequent steps: first, a-particles are formed 
from the free nucleons. Then, in the following a-process, nuclei up to about 
A ps 100 are produced |]. Finally, the neutrino bubble is also an ideal site 
for the r-process synthesizing the elements of about A > 100 jl], ^ [7|. 

The bottleneck of the above nucleosynthesis processes is the formation 
of nuclei with mass numbers A > 9 from the nucleons and a-particles, 
because two-step processes have to be involved in bridging the instability 
gaps at A = 5 and A = 8 (see Fig. 1). Due to the lack of neutrons in 
hydrostatic helium burning in red giants, these instability gaps have to be 
bridged by the triple-alpha reaction 4 He(2a,7) 12 C which proceeds as a two- 
step process a + a <-> 8 Be and 8 Be(a,7) 12 C. However, in the a-process where 
neutrons are available the instability gaps can also be overcome by the re- 
action 4 He(an,7) 9 Be which is much more efficient in the a-process than the 
triple-alpha reaction rate [p], [|. 

Another reaction chain could also be of importance. This chain starts 
with 4 He(2n,7) 6 He and is then continued by 6 He(a,n) 9 Be, thereby closing 
the instability gaps at A = 5 and A = 8 (See. Here one should take into 
account that the abundance of free neutrons for the a-process is of the order 
of 10 20-30 neutrons cm -3 0, |(| and that the yield of the a2n reaction forming 
6 He is proportional to Y a Y£ while the yield of the 2cm reaction is proportional 
to Y^Y n . (The notation Y a and Y n stands for the concentrations of alphas and 
neutrons, respectively.) The rate of the 4 He(2n,7) 6 He reaction was estimated 
in Ref. || and found to be very small. However, as it is shown below the main 
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Figure 1: Possible two-step processes involved in bridging the instability 
gaps at A = 5 and A = 8. The main reaction flow under typical conditions 
in the a-process is indicated by thick arrows, the thinner arrows represent 
almost negligible reaction flows. 



mechanism of the reaction was missed in that study. We thus investigate 
whether the rate of the reaction is enhanced to a sufficient degree due to 
this main mechanism under representative conditions. While we are using 
a three-body model of the process the same issue was studied in a recent 
independent paper |§ using a simpler direct-capture model. 

The 4 He(2n,7) 6 He-reaction is also of special interest because experimen- 
tal data of the nucleus 6 He show features characteristic for halo nuclei ||T0| , [TT|1 . 
This is indicated by the special character of the a-spectra for the decay of 
the J 71 = 2 + level at E = 1.8 MeV in 6 He and the measurement of an abnor- 
mally large electromagnetic-dissociation cross section |T3, ITS]. We note as 



well that the strength functions for the reverse reaction were calculated in 
Ref. jTJJ] in the framework of an approach treating all the three final-state 
particles on equal footing. 



2 The Model 

Let us formulate our model for calculating the 4 He(2n,7) 6 He reaction rate. 
We adopt a two-step mechanism for the process. In the first step the p-wave 
3/2~ resonant state of 5 He is formed and the reaction proceeds through a wing 
of this resonance. Calculations of two-step reactions proceeding through the 
low-energy wing of a resonance are well known from the triple-alpha process. 
However, the mechanism of the reaction we are considering proves to be 
different. In analogy with the triple-alpha process, a transition of a p-wave 
neutron to the resonant state of 6 He followed by a E2 capture to the ground 
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state of 6 He would be the second step of the reaction. However there exists 
an alternative possibility, namely a non-resonant El-capture of an s-wave 
neutron. The latter process proves to dominate strongly the reaction in our 
case. This occurs because of both the lower multipolarity and the absence of 
a Coulomb or a centrifugal barrier suppressing this transition. We thus adopt 
this process as the second step of the reaction. The resonant E2-contribution 
will be commented on below. We use the three-body a+n+n representation 
of the nuclear system in our calculations. 

The reaction rate per particle triplet for 4 He(2n,7) 6 He, in analogy to the 
triple-alpha process, is given by [TBI, 113 



where E\ and E 2 are relative energies in the ( 4 He+n)- and ( 5 He+n)-system, 
respectively. The quantity T ( 5 He, Ei) is an energy-dependent width of 5 He, 
whereas the integrands d (n 4 He) (Ei) /dEi and d (n 5 He) (Ei, E 2 ) /dE 2 are 
transport cross sections a(E)v |l7j for the first and second step, respectively, 
times the corresponding Maxwell distributions ~ \[E exp(— E/kT). 

The first step of the reaction is represented by the first integral in Eq. ([!]) 
[ P~5| , [TEfl . An equilibrium between the production and the decay for the pop- 
ulation of 5 He is assumed here. The relative production rate of 5 He is given 
by ai{Ei)v\. The total production cross section <j\ coincides with the elastic 
( 4 He+n) cross section as neutron emission of 5 He is the only possible exit 
channel. The relative decay rate is given by T ( 5 He, Ei) /h. The integrand is 
the ratio of the 5 He particle density at given E\ to the product of net a and 
n particle densities at the equilibrium. To calculate the energy dependent 
width T ( 5 He, E\) we need a ( 4 He+n)-potential. It was obtained by scaling 
the depth and range of the folding potential to the experimental energy and 
width of the 5 He ground state. Both these quantities are exactly reproduced 
by our potential and cti(£'i) is also reproduced reasonably well. The values 
0.89 MeV [|18| and 0.76 MeV were adopted as the Breit-Wigner 5 He energy 
and width, respectively. 

The value of the width was deduced |19| from the most recent experi- 



ment ||20|| . In Ref. itself a significantly larger value of the width was 
reported and we would like to explain the difference between the two values. 
In Ref. [[2(J the spectrum of that work was fitted by the expression of the 
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form 

^17^72 (2) 

[E na — E Tes — A(E na )] 2 + 



r R (E na ) 



times a smooth function in the (n — a) relative energy E na . Here Tr(E) = 
2-f 2 P(E) and A(E) = --f 2 [S(E) - S(E VCS )} where 7 2 is the reduced width, 
and P and S are the penetrability and the shift function defined as usual. 
The value of Tn(E res ) was then reported as a width. However, this i?-matrix 
value is known to differ from the physical width determining the lifetime of 
a system. Just the latter width T determines the energy-dependent width in 
Eq. (|^) and it can be obtained from the location of the pole E na = E — iT/2 
of the factor (|2|). This width was found in Ref.[19| and the value obtained is 
close to the older estimates |L8|. 



To calculate the second step reaction rate (n 5 He), the 5 He(n,7) 6 He cross 
section is required. As explained above, we consider the electric dipole tran- 
sition El from an incident s-wave. The cross section is given by 

with 

j poo roc 

I=— dn dr 2 u h (r 1 )r l u sc (E 2) r 2 )rlf t=ld=?j/2 (r 1 ,r 2 ) , (4) 
k 2 Jo Jo 

where E 1 = Ei + E 2 + Q\ 2 is the energy of the photon, with Q\ 2 = 0.98 MeV 
being the Q-value of the reaction 4 He(2n,7) 6 He. In Eq. (|]) the coordinates 
7*1 and r 2 are the distances between the valence neutron and the a-core or the 
center of mass of 5 He, respectively. The quantities v 2 and k 2 are the relative 
velocity and wave number in the entrance channel ( 5 He+n) of the second 
step. The quantity Mb( r i)/ r i represents the radial part of the quasi-bound 
wave function of the 5 He resonance, and u sc (E 2 ,r 2 )/r 2 is the radial part of 
the ( 5 He+n) scattering wave function. The function ^=1^=3/2(^1,^2) is the 
radial component in the expansion 



* ( 6 He 



flj( r l> r 2) bUj ('I, <Tlz) Xij (h, V2z)]j =0 ( 5 ) 



of the a+n+n ground-state wave function. Here x&j are angular-momen- 
tum functions with orbital momenta I (necessarily equal to each other) and 
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total angular momenta j = £ ± 1/2, <j\ z and 02 Z are spin variables of the 
valence neutrons, and the brackets [. . .] denote vector coupling. Only the 
(£ = l,j = 3/2)-component may contribute to the cross section. The 6 He 
wave function is normalized to unity and the scattering wave function is 
normalized as w sc — > sin(Ax 2 + S) in the asymptotic region. 

The ground-state wave function Wb^i) of 5 He was obtained as the Siegert 
solution to the ( 4 He+n) Schrodinger equation and this solution was taken in- 
side the centrifugal barrier and normalized there to unity (see e.g. |[21|| ). The 



Siegert solution is an eigenfunction corresponding to the boundary condition 
for a decaying system and inside the barrier it represents the contribution of 
a resonance to continuum wave functions. Its absolute value decreases inside 
the barrier and increases in the outer region. The cut-off radius 7.425 fm 
was taken at the point of a minimum of | Mb | . Replacement of this complex- 
energy (E —iT/2) solution by the real-energy E scattering solution with the 
same cut-off radius leads to less than only 1% difference in the second-step 
cross section. This shows both that the resonance contribution is a predom- 
inant one and that the cut-off radius has been chosen correctly. Variation 
of this radius, e.g., by ±1.5 fm changes the cross section by not more than 
20%. (Reasonable variations of the cut-off radii give even narrower limits.) 
The above-mentioned ( 4 He+n) folding potential was used to obtain u\,. 

To calculate the scattering wave function u sc (E2, r^) we construct the opti- 
cal ( 5 He+n)-potential by using the folding procedure [fj^, ^3j. We first repro- 



duced the matter density distribution for 4 He by means of a simple Gaussian 
with the r.m.s. matter radius of 1.469 fm deduced from the r.m.s. charge ra- 
dius of 1.676 fm [p4j] . To obtain the density distribution for 5 He, we folded the 
4 He-density distribution with the above-mentioned wave function Ub(Ei, ri). 
Finally, we used again the folding procedure to obtain the potential. Due to 
abnormally large separations between the valence neutrons and the a core 
in the final-state 6 He halo nucleus the cross section is largely enhanced. The 
outer free-motion parts of the and u sc functions contribute predominantly 
to Eq. (|4]) due to the same reason. If one replaces the ( 4 He+n) potential 
we used above by the square-well potential whose range and depth are fitted 
to the position and width of the 5 He resonance the reaction rate obtained 
changes by less than 1% and thus it is insensitive to the form of the poten- 
tial. An accurate three-body wave function of 6 He (see p| ) was used in our 
calculations. Originally, it was given in LS'-coupling and in the form of an 
expansion over hyperspherical harmonics and we present it in the form of 
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Figure 2: Calculated cross section for 5 He(n,7) 6 He 



Figure 3: Calculated cross sections for 6 He(7,n) 5 He 

Eq. (|). 

While we included the formation of the 5 He(3/2 ) resonance as the first 
step of the reaction we neglected the contribution from the broad l/2~ res- 
onance. Taking into account Eq. (|l|) the relative l/2~ contribution may be 
roughly estimated as 

[a 1/ 2(E)T 3/ 2Wi/2][ay 2 (E)T 1 /2W 3/ 2)' 1 (6) 

where cr's are elastic cross sections, T's are widths and w's are weights of 
the ft=ij components with j = 1/2 and 3/2 in the 6 He wave function from 
Eq. ([5]). The weights we calculated are w 3 / 2 = 0.87 and Wi/ 2 = 0.049. The 
component entering Eq. (§) thus dominates the 6 He wave function. Sub- 
stituting these values into (||) along with the value of Yi/ 2 ~ 4 MeV ||18 



we obtain 0. Ola i /2(E) / a 3 /2(E). The ratio of the cross sections is less than 
unity. Thus the contribution of the 5 He(l/2~) resonance to the reaction rate 
is negligible. 
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Figure 4: Ratio of calculated reaction rate of 4 He(2n,7) 6 He for Tg = 0.5 — 8.0 
and the maximum of the previous adopted value M 

Table 1: Parameters of the reaction rate for 4 He(2n,7) 6 He calculated in our 
three-body model 





0.1 < T 9 < 2 


2 < Tg < 15 


a 


0.00265 


0.293 


b 


2.55 


-0.351 


c 


0.181 


-5.24 



3 Reaction Rates 

The calculated cross section for 5 He(n,7) 6 He can be parametrized by a = 

1 /2 

0.152/ (E 2 [MeV]) [fib] and it is depicted in Fig. 2. Furthermore, in Fig. 3 
the cross section for the inverse process 6 He(7,n) 5 He, as obtained from the 
principle of time-reversal invariance EH], is shown. The cross sections 
given in these figures were obtained by assuming that 5 He exists exactly at 
the resonance energy E\ = 0.89 MeV. The 5 He(n,7) 6 He cross section exhibits 
the 1/f-behavior corresponding to incident s-waves. The 6 He(7,n) 5 He cross 
section shown in Fig. 3 is about three orders of magnitude higher than the 
maximum value that was originally used for the determination of the reaction 
rates of 4 He(2n,7) 6 He and 6 He(7,2n) 4 He by Fowler et al. [|] (cf. below). 
The obtained reaction rate can be parametrized in the following way: 

Nl (2n 4 He) = a(T 9 ) b exp (Jp) 10" 8 cm 6 s" 1 mol -2 , (7) 

where the parameters a, b, c are listed in Table 1 for the two temperature 
regions 0.1 < T 9 < 2 and 2 < T 9 < 15 (T 9 : in units of 10 9 K). 

The inverse reaction rate A 7 per nucleus per second can be calculated by 
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using the RevRatio 



A 7 = RevRatio x JV£ (2n 4 He) s 



-i 



(8) 



The RevRatio is parametrized in the following way 



RevRatio = 1.08(T 9 ) 3 exp 



( 



-11.3 

T 9 



) 



10 20 cm" 6 mol 2 



(9) 



As shown in Fig. 4, the reaction rate for 4 He(2n,7) 6 He calculated with 
the help of our three-body model is more than three orders of magnitude 
larger than the maximum of the previously adopted value || . This is mainly 
due to the non-resonant El-transition of the second step 5 He(n,7) 6 He, which 
dominates the cross section, and which was not taken into account in Ref. [§]. 
We have also obtained an estimate comparable to the results in Fig. 4 using 
the direct-capture model. 

The huge enhancement of the reaction rate obtained can simply be demon- 
strated as follows. Let us estimate the ratio of the second-step non-resonant 
transport cross section (av) nr to the resonant one. To simplify the reasoning, 
let us consider the case of higher temperatures such that kT is comparable 
to or larger than the width Ti of 5 He (T larger than about 10 10 K). Be- 
low it will be convenient to treat the (crv) nr value as an average transition 
probability. This implies that the 5 He+n relative wave function is normal- 
ized to exp(zkr) plus the scattered wave at large distances. The resonant 
second-step transition disregarded in our calculations proceeds through the 
first excited state 6 He(2 + ) lying within the width of the broader 5 He reso- 
nance. The quantity (av) for this transition can be estimated as (c.f. fl7|l ) 
(^ 2 /2/i) / u~ 1 / 2 A _3//2 a; 3 / 2 exp(— x)hw r where A is the position of the 6 He reso- 
nance with respect to the ( 5 He+n) threshold, x = A/kT, \i is the reduced 
mass for ( 5 He+n) and w r is the transition probability pertaining to the tran- 
sition from the 6 He(2 + ) quasi-bound state to the ground state. The quantity 
A equals to (E r — E\) where E r is the absolute position of the 6 He(2 + ) res- 
onance and E\ is the energy of 5 He. This quantity should be averaged over 
the distribution of E\ giving A(T) ~ Y\. The ratio (<rv) nr /w r of the two 
transition probabilities equals to the squared ratio of absolute values of elec- 
tromagnetic transition matrix elements. The transitions considered are El 
and E2, respectively, and this ratio can thus be estimated as _R 3 [(2/+l)//c 7 i?] 2 
where I = 2, fc 7 ~ 10~ 2 fm _1 is the photon momentum, R ~ 4.5 fm |TT] is the 
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average distance between a ( 4 He+n) subsystem and the other outer neutron 
in the 6 He nucleus and the R 3 factor comes from the difference in normal- 
izations of the initial state wave functions entering the transition matrix 
elements. As a result we arrive at the rough estimate 

-.2 \ -1/2 



25 



\2 



2/i 



Tr 



pR 2 



r 3/ V 3/2 exp(x) ~ 5 ■ 10 3 x" 3/2 exp(x) 



for the ratio of the reaction rates, where x ~ T\/kT. 

The rate obtained in our three-body calculation is about a factor of 1.6 
smaller than the simpler single-particle direct-capture rate given in Ref. M. 
Therefore, we reach the same conclusion as in Ref. ||, i.e. for the trajectories 
relevant to the r-process given in Ref. |4J (the last 16 trajectories given in 
Table III) the reaction rate of the 4 He(2n,7) 6 He obtained is stronger by a 
factor of three than the triple-a reaction, but is still weaker by approximately 
two orders of magnitude than the 4 He(an,7) 9 Be reaction. 

With respect to the photodisintegration of 6 He our calculations show 
that for the same trajectories as given above, the beta decay of 6 He and 
6 He(2n,7) 8 He are slower by about 10 and 11 orders of magnitude, respec- 
tively (the reaction rate for 6 He(2n,7) 8 He is given in [2?|). The reaction 
6 He(a,7) 10 Be is smaller by about two orders of magnitude. The only reaction 
which can compete with the photodisintegration is the reaction 6 He(a,n) 9 Be. 
The cross section for the inverse reaction 9 Be(n,a) 6 He was measured by Stel- 
son and Campbell |28[] in the energy range from 0.7 to 4.4 MeV. From this 
the reaction rate of 6 He(a,n) 9 Be was deduced. A comparison of the rate with 
the photodisintegration rate is given in Table 2 for different temperatures. 
While the rate for the photodisintegration is multiplied with the density p 
in units gem -3 and the 6 He abundance in the network, the rate for the a 
capture is multiplied with the square of the density and the abundances of 
6 He and a. Therefore, under representative conditions (cf., Ref. [l|]) with 
p = 5 x 10 4 gcm -3 and Y a = 0.225 most of the produced 6 He can be pro- 
cessed further on to 9 Be. Even at T 9 = 5.0 alpha capture by 6 He is still more 
important than photodisintegration. However, it has to be noted that un- 
der representative conditions in the a-process the reaction path proceeding 
through 6 He is still negligible compared to the reaction 4 He(an,7) 9 Be. 

In conclusion, we have considered a new mechanism of the two-neutron 
capture by the 4 He nucleus which leads to a more than three order of magni- 
tude enhancement with respect to the value adopted by Fowler et al. B. We 
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Table 2: Comparison of the reaction rates of the photodisintegration 
6 He(7,2n) 4 He (second column) and 6 He(a,n) 9 Be (third column) for differ- 
ent temperatures (first column) 



T 9 


A-y 


N A < 6 Hea > 


0.5 


1.28 x 10" 2 


6.45 x 10 5 


0.8 


7.71 x 10 2 


3.16 x 10 6 


1.0 


4.36 x 10 4 


6.69 x 10 6 


1.5 


1.66 x 10 7 


2.64 x 10 7 


2.0 


4.87 x 10 8 


6.86 x 10 7 


2.5 


4.54 x 10 9 


1.37 x 10 8 


3.0 


2.27 x 10 10 


2.26 x 10 8 


5.0 


8.50 x 10 11 


6.05 x 10 8 



have performed a simplified three-body calculation of the reaction rate that 
proved to be in accordance with the prediction of the simple direct -capture 
model within a factor of two. The inverse reaction of the photodisintegration 
of 6 He has also been studied. Our cross-section for the photodisintegration 
via formation of the intermediate 5 He state can be verified in experiments on 
6 He break-up on heavy targets at low energies, of the type performed last 
years J29], |30| for the n Li nucleus. The results obtained have a relevance to 
other 2N radiative capture processes as well. 
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J.S. Vaagen and M.V. Zhukov for valuable discussions on halo nuclei. Useful 
discussions with J. Gorres, M. Wiescher and M.J. Balbes are acknowledged. 
This work was supported by the Fonds zur Forderung der wissenschaft lichen 
Forschung in Osterreich (project P10361-PHY), by the Osterreichische Na- 
tionalbank (project 5054) and by the International Science Foundation and 
Russian Government (grant J4M100). 
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